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Abstract. According to theory, high energy emission from the coronae of cool stars 
can severely erode the atmosphere of orbiting planets. To test the long term effects of 
the erosion we study a large sample of planet-hosting stars observed in X-rays. The 
results reveal that massive planets (M p sin i > 1.5 Mj) may survive only if exposed to 
low accumulated coronal radiation. The planet HD 209458 b might have lost more than 
1 Mj already, and other cases, like r Boo b, could be losing mass at a rate of 3.4 M e /Gyr. 
The strongest erosive effects would take place during the first stages of the stellar life, 
when the faster rotation generates more energetic coronal radiation. The planets with 
higher density seem to resist better the radiation effects, as foreseen by models. Current 
models need to be improved to explain the observed distribution of planetary masses 
with the coronal radiation received. 



1. Introduction 

The discovery of exoplanets has reached a point in which it is possible to study more 
detailed properties of planets, such as mass evolution. Once the planet is formed, and 
the original disc is dissipated, the main agent interacting with the atmosphere should be 
the high energy emission from the corona of the star, for late type stars. The photons 
with A < 912 A can ionize hydrogen atoms, assumed to be the main component of the 
atmosphere of giant planets. The effects of X-rays (AA 1-100) and Extreme Ultraviolet 
(EUV, AA 100-912) photons take place at different heights in the atmosphere of the 
planet. While EUV photons mainly ionize the atoms in the upper atmosphere, the X- 
rays penetrate deeper in the atmosphere. The free electrons pro duce a cascade of colli- 
sions while the X-rays photons are absorbed in the atmosphere dCecchi-P estellin i et al. 
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2009h . These collisions heat the atmosphere yielding its "inflation" and eventually the 



evaporation of a part of it. The gravity of the planet act as a protection, trying to keep the 
atmosphere attached to the planet. If we assume that the planet atmosphere is mainly 
composed by hydrogen, and all the photons in the whole XUV range (X-rays+EUV) 
are absorbed and contribute to the heating of the atmosphere, it is possible to calcu- 
late the mass loss of the atmosphere by balancing the losses with the planet gravity 
(Wats on et al.lll98 H : lLammer et all 12003). An add itional source of m ass losses takes 



place through the Roche Lobe for close-in planets ( Erkaev et al.ll2007h included in the 



formula through the variable K (K < 1). Some authors consider that evaporation ac- 
tually takes place at a point so mewhere above the pla net radius R p , at the "expa nsion 
radius" R\ dBaraffe et al.ll20040 . The resulting formula (ISanz-Forcada et al1l2010bl) is: 



M = 2 (1) 

GKM p 

where /? - RJR p , Fxuv 1S the X-ray and EUV flux at the planet orbit, and G is the 
gravitational constant. Although some heat may be absorbed at R\, most XUV flux 
should be absorbed below R p , where most of the atmosphere is enclosed. Thus we can 
safely assume ft - 1. The formula can be simplified using the mean density of the 
planet (p), and assume that K - 1 (valid for most cases): 

M=H^V. (2) 
Gp 

Stars with spectral types ranging from mid-F to early-M have stellar coronae with 
temperatures of ~1 MK, exceeding 10 MK in the most active cases. The high tempera- 
ture material in the transition region (~ log T = 4-5.8) and corona (~ log T = 5.8-7.4) 
emits copious X-rays and EUV flux. Fast rotators have a hotter corona, resulting in a 
higher XUV flux. Since the younger stars have a faster rotation, the XUV emission of 
the late type stars decreases w ith time. The evolut i on of X-rays emission with ag e has 
been studied for the Sun (e. g. iMaggio et alJll987l : lAvres!ll997l: iRibas et all 120050 and 



extended to G and M stars (IPenz & Mic ela 2008;: iPenz et al.l l2008). A relation using 
late F to early M stars has been calibrated by Garces et al. (in preparation) in the X- 
rays band, allowing us to calculate the age of the stars from its X-rays emission, and to 
trace the time evolution of this emission. Current instruments allow only access to the 
X-rays band, but there is little information on the EUV emission in stars other than the 
Sun, limited to AA 100-400 for the best cases due to the absorption of the radiation in 
the interstellar medium. 

The distribution of planet masses with current emission from the coronae of the 
stars should reflect the accumulated effects of the mass loss in the atmospheres of the 
planet over time. Besides, if we know the evolution of the emission in the whole XUV 
band we should be able to trace the planet evolution, given an accurate knowledge of the 
density of the planet, according to Eq. |2] We carry out a program to observe the X-rays 
emission of stars with exoplanets and evaluate these and other effects. We have devel- 
oped a database, "X-exoplanets'Q, that includes not only observatio ns of the stars in X 



rays, b ut also synthetic spectra in the whole XUV (1-920 A) range (ISanz-Forcada et al. 



l2010al) . available through the Spanish Virtual Observatory (SVO). We analyzed X-rays 



1 "X-exoplanets" can be accessed from http://sdc.cab.inta-csic.es/xexoplanets 
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Figure 1 . Distribution of planetary masses (M p sin z) with X-ray flux at the planet 
orbit. Filled symbols (squares for subgiants, circles for dwarfs) are XMM-Newton 
and Chandra data. Arrows indicate upper limits. Open symbols are ROSAT data 
without error bars. Diamonds represent Jupiter, Saturn, and the Earth. The dashed 
line marks the "erosion lin e". Dotted lines indicate th e X-ray flux of the younger 
Sun at 1 a.u. Adapted from lSanz-Forcada et al.l d2010bl) . 



data awarded to us either as P.I. or co-I., including an XMM-Newton Large Observa- 
tion (proposals #02065304, #02000 001, #05510229), and arch ival XMM-Newton and 
Chandra data. The scientific results (iSanz-Forcada et al.l l2010b. Sanz-Forcada et al., in 
preparation) are outlined in next section. 



2. Results 

Direct XUV observational data of stars with exoplanets can be acquired only in the 
X-rays band with present instruments. Since X-rays and EUV are formed in the same 
region and temperatures, the trends observed in X-rays should be extrapolative to the 
whole XUV range. We have calculated the X-rays flux by fitting the spectra of all 
the stars with exoplanets observed with XMM-Newton or Chandra, and complemented 
these data with ROSAT observations with S/N>3. A total of 101 planets in 88 stars 
are included in the sample. We can interpret the results keeping in mind Eq. [2] The 



Sanz-Forcada et al. 



23 



FT 



o 

CD 

o 
cE 

Cl 

cb 

(D 



22 



21 



20 



19 



o 



□ 



— o. 



O 



. ■ o u 

C|j» - 

* ® O 



• age > 1 Gyr 

• 10 > age > 1 Gyr 

• 1 > age > 0.5 Gyr 
age < 0.5 Gyr 



o 



CO 

o 



18 



171 



4 6 
Planet mass, Mp sin i (Mj) 



Figure 2. Distribution of planetary masses (M p sin z) with the X-ray flux accumu- 
lated at the planet orbit since an age of 20 Myr to the present (see text). Symbols as 
in Fig.Q] 



distribution of Fxuv with masses (Fig. []} reveal a clear separation that seems to be 
related to masfl We plotted a line ("erosion line") that roughly follows this separation. 
This line is not based on any previous assumption or physical law. The Solar System 
planets are included for comparison, with vertical segments indicating the variations 
over the solar cycle. We also include the radiation level arriving at the Earth when the 
life appeared, 3 .5 Gya, and when the Sun had an age of only 1 Gyr. We used the solar 
young analogs dRibas et al.ll2005l) k 1 Cet and EK Dra to mimic the flux at the Earth's 
orbit in the past. Fig. Q] can be interpreted in the same manner as the HR diagram: the 
absence of planets in a given area of the diagram implies that they spend little time in 
that stage. In this sense the distribution seems to indicate that the planets suffer heavy 
erosion in the first stages until they are below the "erosion line", when either the stellar 
XUV flux decreases and/or the gravity of the planet protects the planet from mass loss, 
slowing down the erosion. Planets positioned initially at further distances will not suffer 
substantial mass loss. We can also divide the diagram in four quadrants based on the 



2 Planets orbiting giants and cataclysmic variables were excluded from the analysis. 
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Figure 3. Density of the 103 planets of known radius (17 Nov 2010) with M p < 
8 Mj (filled circles). Diamonds represent the Solar System planets. 



mass (at 1.5 Mj) and the X-ray flux at the planet orbit (at log Fx=2.15). Only 1 out 
of 12 of the planets receiving high flux have high mass, while 46 out of 89 (a 52%) of 
the planets receiving lower fluxes have a high mass. Therefore it is clear the absence of 
high mass planets suffering high flux levels. 

To test the effects over time we have also used the dependence of X-rays emission 
with time. We need also the EU V evolution with ti me, but the only estimates made 
to date are valid just for the Sun ( Ribas et ai]|2005l) . and a relation applicable to all 
late type stars is necessary. The first step needed is to calculate the EUV flux. The 
process is rather complex and it is better explained in ISanz-Forcada et al". 1 fcOlOah and 
Sanz-Forcada et al. (in preparation). We use the spectral information in X-rays to 
derive a coronal model, extrapolate the model to the transition region if no UV lines 
are available^, and predict the spectral energy d istribution in the whole EUV range 
using the atomic model APED (ISmith et alj|2001l) . We have derived a relation between 
X-rays and EUV luminosity (Leuv ~ L^ 65 ) that we use to calculate the EUV flux of 
the stars with ROSAT data. We also calculated the age of the stars using their X-rays 
luminosity, and constructed a relation in the EUV range with the best data of our sample 
(log Leuv - 29.22-1.41 log t, where t is the age in Gyr). Fig.|2]shows the accumulated 



3 The model of the temperature structure in the transition region is typically calculated from UV lines, 
formed at lower temperatures than those in X-rays. 
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Figure 4. Planetary mass evolution as an effect of the XUV radiation. Mass losses 
include evaporation by coronal radiation and losses through the Roche Lobe. Planet 
mean density (in gem -3 ) is indicated, as well as the current stellar age of the planet, 
calculated using the X-ray luminosity (see text). 



flux at the orbit of the planet, considering the time between an age of 20 Myr and the 
present. The observed distribution with mass would be a direct indication of the mass 
lost by the planets since an age of 20 Myr, assuming that all mass losses are included in 
Eq. |2j and the density is the same for all planets. The distribution shows that only three 
planets with M p > 1.5Mj have survived an X UV flux of more than 10 21 erg cm -2 : for 
two of them there are alternative explanations dSanz-Forcada et aU2010b ) and the third, 
t Boo b, is a young object that might still be suffering heavy erosion. Alternatively, if 
any of these planets has a rocky composition (p >5gcm~ 3 ) the erosion effects should 
be much lower, if any. 

The second important variable in Eq. [2] is the planet mean density. A more dense 
planet is better protected against evaporation. Ideally our analysis should concentrate 
on stars with known Lx and planets with known density. Only four planets have both 
magnitudes known, and one of them, 2M1207 b is too far from the star to suffer mass 
losses. Instead we can test the whole population of planets with known density, a total 
of 103 planets. Most of these planets have been detected trough the transits technique, 
therefore having a strong bias towards close-in planets, those more exposed to high 
XUV flux. We do not find massive planets with small density among this population 
(Fig. |U). This might indicate that low density planets with high mass have suffered 
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a quick erosion taking them to the population with jovian masses, or increasing their 
density because of the lower atmosphere-to-nucleus ratio. 

Finally we have made a calculation of the mass evolution for the mentioned three 
planets with known density and Lx, including the calculated Levy, according to Eq. |2j 
and including in this case the mass losses through the Roche Lobe expressed in Eq. Q] 
Fig. |4] shows this evolution assuming that density has been roughly constant over time 
and no other effects toke place. The case of HD 209458 b is remarkable, with more 
than 1 Mj lost to date. This work has considered only a rather simple model with 
the thermal losses in the atmosphere, but non-thermal losses should be included in 
the future, together with other smaller effect s that may be nece s sary to explain the 
distribution observed in Fig. [T] as described in lSanz-Forcada etaL (I2010bh . 
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